Structural Stability of Vacancy-ordered YD2.75C60 
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The fcc-based structure of Yb2.7sC6o is unique among metal-doped fullerene compounds, exhibiting 
long-range-ordered vacancies, significantly off-centered divalent Yb cations, and distorted, crystallo- 
graphically inequivalent, orientationally ordered Ceo anions. A simple electrostatic-energy analysis, 
which models the constituents using point charges, is shown to provide insight into how each of 
these features stabilizes this unusual crystal structure. The results have general implications for a 
variety of other intercalated metal fullerides. 

PACS numbers: 61.48. +c, 61.50.Ah, 74.70.Wz 

Since 1991, numerous superconducting compounds of 
Ceo intercalated with monovalent alkali and divalent 
alkaline-earth metals have been prepared and struc- 
turally characterized jjj. In many respects, their crystal 
structures can be described by regarding the Cgo anions 
as large, negatively-charged "atoms" forming simple fee 
or bec lattices, the interstices of which are filled by metal 
cations. The stability of these structures is generally well 
understood within simple models dominated by electro- 
static energies [|],[| . The analysis of the structural ener- 
getics through more complete first-principles calculations 
confirms and refines the simple models Q . 

Ytterbium metal has also been shown to form a su- 
perconducting compound with Gqq, but with a structure 
that is considerably more complex jH) . It can be viewed as 
a modification of a hypothetical YbsCgo structure, analo- 
gous to that of K3C60 H and schematically illustrated in 
Fig. [ija). In this ideal Yb 3 Ceo structure, Yb cations sit 
at the centers of the octahedral (O) and tetrahedral (T) 
holes created by the fee lattice of pure Ceo- Systemati- 
cally removing 1 of the 8 T-cations leads to the ordered- 
vacancy structure of Yb2.7sC6o in Fig. |l](b), where the 
full unit cell is seen to contain 8 subcells, each with a 
different T-site vacancy. One such subcell, with the Ceo 

FIG. 1. (a) Schematic picture of hypothetical fee YbaCeo 
structure (Ceo anion radius is 3.5 A, Yb 2+ radius is 1.0 A). 
Cations fill all tetrahedral (T) and octahedral (O) interstitial 
sites and are shaded gray and black, respectively. Removing 
one T-site cation (shown with arrow) leads to the fcc-based 
vacancy-ordered Yb2.7sC6o structure in (b), with the Ceo an- 
ions removed for clarity. Eight subcells comprise the full unit 
cell. In the expanded subcell shown, O-site cations are signifi- 
cantly displaced from their ideal, "centered" positions towards 
the nearest vacancy, and each T-site cation is displaced less 
(drawn exaggerated) along one of the Cartesian axes. Other 
O-site cations have been displaced towards vacancies in adja- 
cent subcells. (c) Same Yb2.7sCeo subcell as in (b) but with 
the Yb cations omitted. The vacancies create 3 inequiva- 
lent types of Ceo anions, schematically shown with different 
shadings and shapes. Each anion rotates about the local axis 
indicated, to maximize the number of pentagon faces oriented 
towards the surrounding cations. 
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anions removed for clarity, highlights the large displace- 
ments of the O-cations towards their nearest-neighbor 
(NN) vacancy and the smaller displacements of the lo- 
cations along Cartesian axes f|. The 3 types of crys- 
tallographically inequivalent Ceo anions with 2, 1, or 
neighboring vacancies, and their correspondingly differ- 
ent coordinations of 12, 10, and 8 Yb cations, are shown 
in Fig |l](c). Also included schematically are the differ- 
ent distortions of the anions and the local axes about 
which they rotate in order to orient their pentagonal faces 
towards the NN cations. 

In this Letter, we identify the underlying reasons for 
the stability of this apparently exotic phase. A simple 
electrostatic model is used to evaluate the energies as- 
sociated with vacancy formation, cation displacements, 
and Ceo orientations and distortions. Our analysis shows 
that the latter relaxation energies more than compen- 
sate the cost of creating a Yb vacancy. The formation of 
Yb2.75C6o is therefore seen to fit naturally into a unified 
understanding of both the stability of metal-doped Ceo 
compounds and the role of particular structural defects. 

Our investigation focuses on understanding why the 
ideal YbaCgo structure in Fig |l|(a) is not formed but 
is instead found to phase separate into the observed 
Yb2.75C6o structure and Yb metal. For these two struc- 
tures, we compare the relative energies of steps in a Born- 
Haber cycle, taken with respect to a reference system of 
undoped fee Geo and isolated neutral Yb atoms. The 
steps are (i) ionizing the Yb atoms into Yb 2+ + 2e~ ||, 

(ii) condensing any excess Yb atoms into Yb metal, (iii) 
completely transfcring the ionized electrons from Yb to 
the Ceo molecules, and (iv) assembling the cations and 
anions into a crystal. Since the electrostatic energy is 
the largest single contribution to the structural energetics 
in these systems, it is expected to dominate the energy 
differences between the configurations of interest. Ac- 
cordingly, we approximate the relative energies in steps 

(iii) and (iv) by the electrostatic energy differences alone. 
The ions are modeled using point charges, and the Ewald 
method || is used to evaluate the resulting Coulombic 
sums. Experimental structural parameters are used to in- 
corporate the effects of short-range repulsive interactions 
between the Yb cations and Cgo anions. The parameters 
from Ref. 5 are slightly simplified by taking a = b = 
c = 27.8733 A (i.e., ignoring the < 0.4% orthorhombic 
distortion) and by assuming complete occupation of the 
octahedral sites. All energies are given in eV per YbaCeo 
formula unit, or equivalently, in eV/Cgo- 

We begin by computing the relative formation energies 
of two hypothetical structures, ideal YbsCgo and "unre- 
laxed" Yb 2 .75C 60 . The latter is obtained by removing 
the requisite Yb T-cations from an ideal YbaCgo lattice 
to form the ordered- vacancy structure in Fig 0(b), while 
keeping all the remaining cations at the centers of their 
respective O- and T-site holes. The ionization steps for 
ideal YbsCeo and unrelaxed Yb2.7sC6o involve, respec- 



tively, 3 Yb atoms with an energy Ei on = 55.29 eV and 
2.75 Yb atoms with E ion = 50.68 eV |Q. The excess 
0.25 Yb atom per Cgo in Yb2.7sC6o is condensed into 
Yb metal, giving a cohesive energy gain Tv co h = —0.4 eV 
pl| . The total electron affinity for transferring 6 elec- 
trons to Ceo in YbsCgo is approximated as 6 times the 
first electron affinity for Cgo jl2| plus the electronic self- 
energy, which is modeled as the minimum energy of 6 
classical point electrons on a uniform spherical shell of 
radius 3.5 A Q; the result is a net energy cost E c& = 
-15.90 + 41.04 = 25.14 eV. The analogous energy for 
transferring 5.5 electrons to Ceo in Yb2.7sC6o is approx- 
imated by interpolating the self-energy between integer 
numbers of electrons n = 2, 3, 4, and 6 with a second- 
order polynomial in n, yielding E ca = —14.58 + 33.17 
= 18.59 eV. Finally, the Yb cations and Ceo anions are 
assembled into their respective fec-based crystal struc- 
tures, differing only by the T-site vacancies in unrelaxed 
Yb2.75C6o- The Madelung energies, computed by treat- 
ing the Yb cations as points with charge 2+ and replac- 
ing the spherical shells of the Cgo anions with points of 
charge 6— or 5.5—, are Emm = —91.43 eV for YbsCgo 
and —77.65 eV for Yb2.75Cgo- Comparing the total for- 
mation energy E fmm = E ion +E coh +E c& +E M M = -11.00 
eV for ideal YbsCeo with Ef olm = —8.78 eV for unre- 
laxed Yb2.75C6o shows that the larger Madelung energy 
gain in forming YbaCeo exceeds its larger ionization and 
self-energy costs, and thus creating T-site vacancies in 
unrelaxed Yb2.75Cgo is endothermic: ATf orm = 2.22 eV. 

The possible relaxation-energy gains associated with 
vacancy formation must now be considered. One obvi- 
ous source is the displacement of each O-site Yb cation 
towards its NN vacancy (see Fig 0(b)). The equilib- 
rium value for this displacement, determined by the bal- 
ance of electrostatic forces and short-range repulsion be- 
tween the cations and their neighboring anions, is known 
from experiment to be extremely large (2.41 A |||); us- 
ing this value leads to a computed energy gain of —2.32 
eV. The much smaller off-center displacements of the T- 
site cations (~0.3 A) contribute just —0.009 eV on their 
own, but combined with those of the O-site cations lead 
to a small cooperative interaction and a total relaxation- 
energy gain T disp i = -2.40 eV. 

The above calculations assumed that the charges on all 
Ceo anions are given by the average number of electrons 
per Cgo, namely — 5.5e (^]. The crystallographic inequiv- 
alence of the anions (see Fig 0(c)), however, makes favor- 
able the prospect of transferring electrons from Cen(l), 
with 12 NN cations, to C 60 (S), which has only 8 @. We 
investigate this possibility by allowing the inequivalent 
Cqo(1), Cqo(2), and Cqo(3) anions to have correspond- 
ingly different charges — q^e, —q2&, and —qge, with the 
only constraint being that {qi+6q2+q3)/8 — 5.5. Mini- 
mizing the quantity i?Mad + T Ga with respect to qi, we 
find qi — 6.3, 172 — 5.5 and qs = 4.7. The charge of Cgo (2) 
is unchanged, while 0.8e~ is indeed transferred from 
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Ceo (3) to Ceo(l)- This is not a charge disproportion- 
ation but is rather a charge redistribution, reminiscent 
of the different charge states of the crystallographically 
inequivalent planar and chain O anions in YBa2Cu307 
|7|,[l4|. We calculate that in Yb 2 .75C6o the associated 
energy gain -E re dist is just —0.17 eV. This small value re- 
flects the nonlinear self-energy of electrons on Ceo, and 
shows that the redistribution process plays only a minor 
role in determining structural stability. The magnitude 
of the relaxation energy gain in Yb2.75C6o of -Edispi + 
-Eredist = —2.57 eV is, within the estimated accuracy of 
our calculations, essentially equal to A_Ef orm = 2.22 eV. 

We now consider the additional relaxation energy that 
can come from optimizing the orientation of the Cgo an- 
ions with respect to their anisotropic environment of NN 
cations. The largest contribution to the orientation en- 
ergy is the electrostatic interaction between the charge 
density of electrons transferred to Cgo, anisotropically 
distributed over its surface, and the neighboring Yb 2+ 
cations jlB). To calculate this orientational interaction 
energy obviously requires going beyond approximating 
the Cgo anions as uniformly charged spherical shells. 
Based on comparisons with chemical analogs [Q,[l5],[l6| 
and first-principles calculations of Cg^ 0, which show 
that the transferred charge density is primarily associ- 
ated with the pentagons, we model a Cgo anion as a 
truncated icosahedron with equal point charges localized 
at the centers of its 12 pentagonal faces and positioned 
3.5 A from its origin |L7|]. For a single cation neighbor, 
the Ceo anion will orient itself to minimize the distance 
between one of the point charges and the cation, mim- 
icking the preference for a cation to be over a pentagonal 
face. Such an attraction between Cgo pentagonal faces 
and neighboring cations was first observed in (Ba 2+ )3CgQ~ 
|l8f and was found to be an important component in 
achieving the best Rietveld refinement of the Yb2.75C6o 
structure IHUlSfl . 

The orientational interaction energy, Eq, was investi- 
gated using the set of {111} axes established in the x- 
ray structural determination || ; these axes are shown in 
Fig. 0(c). Each C 60 anion is rotated by the same angle 
around its own local axis. O is taken to be zero in the 
uniform reference configuration where the twofold axes of 
the Cgo anion are aligned along the Cartesian axes [po[ . 
As shown in Fig. 2, E® for both Yb 3 Ceo and Yb2.7sC6o 
depends strongly on O and exhibits a pronounced mini- 
mum at O = 37.5°, consistent with the x-ray refinement 
results ||[l9]]. However, because of the large displace- 
ments of the O-site cation in Yb2.75Cgo, a greater degree 
of pentagon-cation interaction is achieved through this 
orientation of the Cgo anions. Relative to ideal YbsCgo, 
the orientational relaxation -energy gain in Yb2.7sCgo is 
appreciable: AE@ — —0.97 eV. While we expect our sim- 
ple point-charge model for distributing charge on Cgo to 
overestimate the absolute magnitude of the orientational 
dependence of the energy, the scale of the calculated dif- 
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FIG. 2. Angular dependence of orientational interaction 
energy in the Yb2.7sC6o and ideal YbaCeo structures due to 
rigid rotation of all Ceo anions about local {111} axes. 

ference AEq demonstrates the importance of this effect 
for stabilizing the observed vacancy-ordered phase. 

A final source of relaxation energy can arise from dis- 
torting the Cgo anions. The anisotropic crystal field 
of NN Yb cations (see Fig|l](c)) is expected to distort 
the symmetric shape of a Cgo anion even with a uni- 
form charge distribution, and such distortions will be still 
greater when the charge distribution is itself anisotropic 
P,[7 yT^JT^ |. Indeed, recent x-ray absorption and Raman 
measurements in Yb 2 .7sCgo show clear evidence for sig- 
nificant shape deformations of the Cgo anions 0. We 
calculate the relative distortion -energy gain A-E-dis using 
the same point-charge model as in the calculation of E<~> , 
with the charges allowed to move radially relative to their 
ideal positions in the = 37.5° orientation. The stiffness 
of the Cgo anion is modeled by connecting the charges 
to the anion center with harmonic springs whose spring 
constant is chosen to reproduce the observed A g (l) radial 
breathing mode frequency of ~500 cm" 1 |2l}| . The associ- 
ated relaxation-energy gain is —0.32 eV in Yb2.7sCgo and 
-0.28 eV in Yb 3 Cg , giving AE dis = 0.04 eV. In both 
systems, outward displacements range from 0.01 A to 
0.04 A. While these displacements are somewhat smaller 
than those inferred from the EXAFS data , suggesting 
that the absolute distortion energies are underestimated 
by our model, the scale of the difference AE^is indicates 
that the distortion contribution to the total relaxation 
energy is relatively minor. 

Our comparison of relevant energies in forming the 
Yb2.75Cgo and YbsCgo structures provides strong sup- 
port for the central role played by the interaction be- 
tween the anisotropically positioned Yb cations and the 
anisotropically distributed charge on the Cgo anions. As 
shown here, this strong interaction can be investigated 
through a relatively straightforward analysis of electro- 
static energies, and does not require invoking any chemi- 
cal behavior peculiar to Yb. Despite its simple form, the 
interaction has a profound effect on compound forma- 
tion and crystal structure. In particular, application of 
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our model gives insight into the stability of a seemingly 
unwieldy structure involving both long-range-ordered va- 
cancies in a 2008-atom unit cell and crystallographically 
inequivalent Ceo anions with low-symmetry cation coor- 
dinations. The relationship between these unusual struc- 
tural features is now clear: A T-site cation vacancy leads 
to the O-site cation displacements, which in turn amplify 
the orientational cation-anion interactions. Only one va- 
cancy per subcell is needed to induce displacements of all 
O-cations, making the formation of additional vacancies 
unfavorable. This naturally accounts for the observed 
Yb2.75C6o stoichiometry. 

One consequence of this simple picture is that the same 
ordered- vacancy phase for Yb2.75Cgo should be observed 
in fulleride compounds containing other divalent cations 
whose ionic radii and elemental cohesive energies are sim- 
ilar to Yb. Such expectations are, in fact, realized for 
both Ca- HH and Sm-doped §§ C 60 . Since the mag- 
nitudes of the electrostatic energies scale quadratically 
with charge, it might also be expected that many of the 
stabilizing factors identified here for these divalent-cation 
fullerides should occur for the similarly-sized monovalent 
Na cation. However, neither a NaaCgo structure nor a 
vacancy-ordered fulleride phase is observed; instead, a 
structure in which two Na cations occupy the same O- 
site hole is found (24). The reason for this is that other 
factors, such as the short-range repulsive interactions be- 
tween cations and Gqq anions and between cations in the 



same O-site 25 1, scale differently with charge and affect 
the structural stability of this system. 

The question of off-center displacements for O-site 
cations has recently been a subject of debate in several 
fcc-based alkali fullerides |2^,|2^]. For those systems in 
which T-site vacancies are observed, e.g., Rb 3 Ceo p8j , 
our analysis leads unavoidably to the conclusion that the 
surrounding O-site cations will undergo off-center dis- 
placements towards that vacancy (see Fig 0(b)). While 
this applies only to the fraction of O-site cations im- 
mediately surrounding the vacancies, such displacements 
should be detectable with appropriate measurements. 

The model presented here is designed to capture the es- 
sential features of the structural energetics of the systems 
of interest, and to enhance the understanding of the ori- 
gin and properties of complex structures. In Yb2.7sC6o, 
we found that the relaxation energies associated with the 
lowered symmetry are crucial in stabilizing its formation. 
The quantitative limitations of this type of analysis are 
clear. However, our model can provide valuable guidance 
for future, more definitive investigations of structural en- 
ergetics in this and similar systems using first-principles 
calculations, and for designing and interpreting related 
experiments in the interim. 
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